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Abstract A new 3D transition turbulence model, more accurate and faster than 
an empirical transition model, is proposed. The model is based on the calculation 
of the pre-transitional u v' due to mean flow shear. The present transition model is 
fully described and verified against eight benchmark test cases. Computations are 
performed for the ERCOFTAC flat-plate T3A, T3C and T3L test cases. Further, 
the model is validated for bypass, cross-flow and separation induced transition and 
compared with empirical transition models. The model presents very good results 
for bypass transition under zero-pressure gradient and with pressure gradient flow 
conditions. Also the model is able to correctly predict separation induced tran¬ 
sition. However, for very low speed and low free-stream turbulence intensity the 
model delays separation induced transition onset. The model also shows very good 
results for transition under complex cross-flow conditions in three-dimensional ge¬ 
ometries. The 3D tested case was the 6:1 prolate-spheroid under three flow condi¬ 
tions. 
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1 Introduction 


This work is the result of a continuous research effort on turbulence transition 

[2013b] 


models development, Vizinho et al. 


Vizinho et al. [2013a) and Vizinho 


et al. [2012]. Transition modeling is still not widely used in industrial flow compu 


tations, mainly due to several limitations of present day models. This seems to be 
the case since the available transition models tend to be complex. Thus, compu¬ 
tation time is largely increased by the usage of a transition tool just to calculate 
an often small fraction of the relevant flow. It goes without saying that it is of 
vital importance to know where transition occurs. However, for some engineering 
problems this transition onset is overlooked and fully turbulent models are applied 
instead, Pascoa et al. [2006 and Pascoa et al. [2010 . Some of these, specifically 


low-Reynolds turbulence models, are sometimes used to predict transition onset. 
The few first serious efforts to evaluate transition modeling capabilities of used 


turbulence models were presented by 

Savill 

[1995 

. Although some of these models 

have an apparent transition behavior, 

Rumsey 

2007 presented a work showing 

that this behavior is mere coincidence and can 1 

ead to design mistakes. It is then 


of great interest to have simple robust transition onset prediction tools. 

The use and general interest in turbulence transition modeling and control in 


industrial engineering is steadily increasing, as seen in the works of Aupoix et al. 
[2011] , Langtry and Menter [2006 , Xisto et al. [2012] , Xisto et al. [2013 and Ab- 
dollahzadeh et al. [2014] . This also applies for the MAAT project presented by 
Ilieva et al. 2012 , Trancossi et al. [2012 and Ilieva et al. [2014] , a revolutionary 
concept of transportation. Due to its high altitude flight conditions, low turbulence 
intensities are predominant. The delay of transition to turbulence is inversely pro¬ 
portional to the free-stream turbulence intensity. This implies that for the MAAT 
project vehicles, feeders and cruiser, large regions of laminar flow will occur. Thus, 
elements of the feeders and cruiser vehicles, such as propulsive systems, will be ex¬ 
posed to late transition and consequently considerable sized laminar flow regions. 
Possession of a turbulence model that can accurately predict transition onset, as 
well as correctly compute the transition length, has become a must amongst CFD 
engineers who deal with transitional flows in a daily basis. 

A common approach used to tackle these flows was the definition of laminar and 
turbulent regions on a priori knowledge of the transition zones. Though functional, 
this approach could only compute turbulent flow over the extension previously 
defined by the user. Also, transition length was not taken into consideration in 
this method. This methodology was impractical in prototype cases, since there is 
no initial information regarding transition onset. 

Then, through the simultaneous work of Smith and Gamberoni [1956 and In 


gen 


1956 


the e n method was developed. Considered by some to be the state of 
the art in prediction of transition onset, it is a method based on the linear local 
stability theory. According to van Ingen [2008], “The linear stability theory consid¬ 
ers a given laminar main flow upon which small disturbances are superimposed”. 
The operation protocol of this method begins with the resolution of the laminar 
boundary layer velocity profiles. This is then used to compute the growth rate of 
the superimposed linear instabilities within the flow by solving the linear stability 
equations. After completion of the latter, an integration of the obtained growth 
rates throughout the streamlines is performed. The resulting value will represent 
the amplification factor, n, of the instabilities. The disturbance amplitude ratio 















































































Transition V-model 


3 


is then given by e". When this variable reaches a threshold, transition onset is 
assumed. Typical transition threshold values for e n are obtained with n values 
between 7 to 9. 

Alternatives for transition prediction were proposed, such as in the experimen¬ 
tal work of Schubauer and Skramstad [1948 that resulted in one of the first em¬ 
pirical graphic correlations for transition onset. Later, the works of |Abu-Ghannam| 
and Shaw [1980] , Mayle [1991] and others produced various empirical correlations 
for transition onset. Based on these empirical correlations, transition to turbu¬ 


lence models were formulated such as the ones created by Cho et al. [1993 , Suzen 


and Huang [2000 and Steelant and Dick [2001] . Typically these first empirically 


correlated transition models were non-locally formulated, since they depended on 
integration of the boundary layer velocity profiles. The latter was performed so 
as to obtain one of the parameters used by the empirical correlations. Most com¬ 
monly this would be the local momentum thickness Reynolds number. Recently, 
according to the work of Menter et al. [2002] and |Langtry| [2006] , locally formu¬ 
lated empirical correlation transition to turbulence models have been developed. 
The previous cited work presented a transition model able to compute transition 
onset as well as transition length. The improvement over the previous empirical 
transition models is that in this formulation the variables are locally defined. In 
Langtry 2006 a transition to turbulence closure is coupled with the SST-k-a> tur¬ 


bulence model of Menter [199 1 . Although, according to the former, it is possible 
to use these transition components of the model with other turbulence models. 
The model formulation can be resumed down to two transport equations for the 
transition part of the model. One for 7 , or intermittency, and another for Re^, 
the transition momentum thickness Reynolds number. A threshold for momentum 
thickness Reynolds number is first computed in the free-stream. This is done ac¬ 
cording to the employed empirical correlation. This value is then diffused into the 
boundary layer. The momentum thickness Reynolds number is computed locally 
with an algebraic relation using a vorticity Reynolds number as presented in the 
works of Van Driest and Blumer [1963] and Menter et al. [2002 . When the lat¬ 
ter computed value reaches the diffused transition threshold momentum thickness 
Reynolds number, transition onset is assumed. The production term in the inter¬ 
mittency transport equation is activated. This in turn is coupled to the selected 
turbulence model. Its production and destruction terms are multiplied by the effec¬ 
tive intermittency allowing control of transition onset. The main problem with the 
work of Langtry [2006] was the lacking of two key components of the model, due to 
proprietary reasons. However and thanks to the work of Suluk sna and Juntasaro| 
[2008] , Suluksna et al. [2009] and Paul Malan and Juntasaro [2009] these lacking 
terms were approximated by alternative functions. Nonetheless, later in [Langtry 
and Menter [2009], the original authors of the 7 — R e e transition model published 


all of the lacking terms. This transition model has a very promising future as an 
engineering tool. Also, the work of Durbin [2012] presents a transition prediction 


tool based on turbulence intermittency that is not empirically correlated. It is an 
intermittency transport equation that does not rely on external data input, such 
as empirical correlations. Similarly to the latter described model, this equation 
can be coupled to any turbulence model making it a versatile tool for transition 
prediction. 

Turbulent flow has always been distinguishable from its laminar state. Transi¬ 
tion can be identified as an increase in skin-friction and velocity profile departure 
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from a Blasius distribution. There use to be a certainty about the fact that only 
turbulent flow had velocity fluctuations. Laminar flow was considered to develop in 
a stationary fashion. This all changed with the work of Dryden 1937 . Considered 
to be the first attempt to successfully measure velocity fluctuations in the laminar 
flow region. It was concluded by Dryden 1937 , that it is impossible to perceive a 


distinction between laminar and turbulent flow solely based upon velocity fluctua¬ 
tion measurements. These fluctuations in the laminar flow region were attributed 
to the presence of a turbulent free-stream. Later on, the experimental work of 
Schubauer and Skramstad [1948 also confirmed the presence of laminar fluctua¬ 


tions just before turbulence transition onset. In the latter experiment, free-stream 
turbulence intensity was reduced as low as possible. Contrary to the previously 
recorded laminar fluctuation values from Dryden 1937 , it was expected that in 
this particular case, the velocity fluctuations would be greatly reduced. Since these 
were related to the free-stream turbulence intensity, the results should bear almost 
no trace of laminar fluctuations. This was confirmed for the leading edge regions of 
the laminar flat-plate flow. Though, the progressive measurement of fluctuations 
along the flow direction began to detect weak oscillations in the laminar region. 
These increased towards the transition up to the bursting of turbulent spots. To 
the best of our knowledge, and according to the award winning paper of |Mayle| 
[1997] , the work of Lin [1957 was the first to analytically evaluate 


and Schulz 


the effects of laminar fluctuations over laminar velocity profiles. This study con¬ 
firmed the possibility of having velocity fluctuations in a laminar flow, even when 
maintaining a Blasius velocity profile distribution. This fact was also documented 
by Dyban et al. [1976] , |Sohn and Reshotkoj [1991] and |Zhou and Wang] [1995] 


In 


the work of Mayle and Schulz [1997] the LKE, or, Laminar-Kinetic-Energy theory 
was proposed. Following this development, a new tool for transition modeling was 
available. The first transition models based on LKE were developed by La rdeau] 
et al. [2004| and [Walters and Leylek| [2004| . Others followed this trend, such as 


Vlahostergios et al. [2009] . In the publication of Lardeau et al. [2004] , it is stated 
that, before the increase of skin-friction due to transition, a growth of fluctua¬ 
tion intensity persists in the upper to medium regions of the laminar boundary 
layer. Although these oscillations have their origin in the free-stream turbulence, 

they do not have the same known fully turbulent ratio of ~ u k v ~ 0.3. Instead 
they present far lower values than the latter. As discussed in |Mayle and Schulz| 
[1997], these fluctuations do not belong to a normal turbulent regime. These are 
then the laminar kinetic energy fluctuations that appear in the upper medium 
regions of the boundary layer. These streamwise fluctuations are Klebanoff modes 
identified by Klebanoff [1971 . The production of these fluctuations have particu¬ 
lar characteristics. The boundary layer selectiveness of free-stream turbulent eddy 
scales filters the broad spectrum free-stream turbulence. This has been identified 
as “shear-sheltering“ and was first described by Hunt et al. [1996 . The work of 
Jacobs and Durbin [1998] , demonstrates this effect using both the discrete and 


continuous modes of the Orr-Sommerfeld equation. The continuous modes are 
eliminated from the main boundary layer region, being confined to the upper 
reaches of the latter. It was also concluded that the boundary layer penetration 
depth by the continuous modes is inversely proportional to their frequency. This 
enforces the concept that only low frequency disturbances are amplified by shear 
in the pre-transitional laminar boundary layer region. In the experimental work 

































































Transition V-model 


5 


of 


Volino and Simon 


1994 


spectra of fluctuating streamwise velocity u , wall- 


normal velocity v and tur bulent shear stresses — u v were recorded. It was found 


in that work and also in 


Volino and Simon 


1997b 


noticeable values of —uv 


in the pre-transitional boundary layer. These were correlated with peak values of 
low-frequency wall-normal velocity fluctuations v of the free-stream turbulence. 
As previously mentioned these — u v values had lower energy and frequency than 
those found in fully turbulent boundary layer flow. As discussed in the work of |Leib| 
Klebanoff [1971] , first named these streamwise velocity fluctuations as 


et al. 


’’breathing modes 11 . Taylor 1939] , noticed that these oscillations were related with 


thickening and thinning of the boundary layer. The production of these stream- 
wise velocity fluctuations u , are believed to be related to the wall-normal veloc ity 
oscillations v through the ” splat-mechanism “ mentioned by 


Bradshaw 


1994 


or 


by the concept of ” inactive mot ion “ p roposed by Townsend [1961] and Bradshaw 


1967 


As explained by 


Volino 


1998 


a negative v velocity fluctuation imposed 


by a turbulent eddy will momentarily compress the boundary layer, shifting higher 
speed flow against the wall surface. This results in an increment of u . As the tur¬ 
bulent eddy is convected by the flow, the imposed compression effect is diminished 
resulting in a recovery of the boundary layer to its previous state. 


The proposed transition model presents some behaviors of the just described 
processes. The most noticeable is the prediction of low negative values of uv' in 
the upper regions of the pre-transitional boundary layer as will be later shown. 
Also the model predicts when these values of uv' pierce the laminar-boundary 
layer. This is then known as the transition onset. 


The main purpose of this work is to present the rational behind the devel¬ 
opment of a new transition model, henceforward designated as V-model. The 
transition V-model is coupled to the Spalart-Allmaras turbulence model and is 
designated throughout the present work as V-SA. First, a short description of the 
three main turbulence transition mechanisms will be presented. The general V- 
SA model coupling is disclosed. Based on the highlighted physics of transition we 
then proposed the new mechanical model equations. The description of the used 
mechanical model approximation is done by first considering some pre-transitional 
turbulent kinetic energy relations. Afterwards the mechanical model approxima¬ 
tion is presented in detail. The development of the transport equation for the 
pre-transitional turbulent kinetic energy of the proposed transition model is pre¬ 
sented. Finally the detailed transition V-model and Spalart-Allmaras turbulence 
model coupling is described. There will also be a validation of the model under var¬ 
ious geometries and turbulent flow conditions. Initially the model is validated for 
zero-pressure-gradient bypass transition using the flat-plate T3A test case. Then 
the model is validated for pressure-gradient bypass transition using the flat-plate 
T3C3 test case. Thereafter the model is validated for separation induced transi¬ 
tion by using the flat-plate T3L test cases. Finally a validation of transition under 
cross-flow effects is performed over the three-dimensional 6:1 prolate-spheroid ge¬ 
ometry. 
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|Pre-Transition| [Transition] iTurbulent" 


Fig. 1 General overview of laminar boundary layer transition to t urbulence. Som e of the 
concepts of late transition depicted here were observed in the work of|Liu and Chen|[2011 . 


2 Transition mechanisms 

Generally speaking we may define diverse transition mechanisms, according to the 
physics of the flow. In this work we develop a model that focuses on the pre- 
transitional region depicted in Figjl] 


2.1 Natural transition 

For flow with free-stream turbulence intensity, FSTI< 1%, natural transition is 
generally observed when the developing laminar flow reaches a critical Reynolds 
number value. After this critical stage of flow development, viscous Tollmien- 
Schlichting instability waves begin to slowly grow, from small linear perturbations 
to non-linear disturbance waves. Having reached this point, the critical location is 
usually defined where the first instability starts to grow, and the transition point 
where the first turbulent spot appears in the transition process. The non-linear 
waves create three dimensional disturbances and, by means of inviscid mecha¬ 
nisms, spots of turbulence begin to appear randomly inside the laminar boundary 
layer. Some of these spots grow in size and downstream from where transition first 
began to develop, a process of fusion between these spots gives light to the fully 
turbulent boundary layer. 


2.2 Bypass transition 

Bypass transition may be sub-critical, it may occur upstream of the previously 
defined critical location. Bypass transition normally occurs due to two sources of 
disturbances, surface roughness and high FSTI, which translates to FSTI> 1%. 
The bypass transition mechanism can be shortly described as a natural transi¬ 
tion process without the Tollmien-Schlichting wave development phase. Therefore, 
laminar flow under bypass transition will turn from laminar to the turbulent spot 
surge phase immediately. The subsecprent development is similar to the natural 
transition. 
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Fig. 2 V-SA model architecture, by coupling the transition V-model to the Spalart-Allmaras 
turbulence model. 


2.3 Separation induced transition 

Separation of flow is generally observable in flow development over surfaces under 
adverse pressure-gradient conditions. Due to flow separation a vortex bubble is 
formed and flow may, or may not, reattach by closing the bubble. The reattach¬ 
ment process is dependent on the increase of flow mixture capacity originated by 
turbulence. Laminar flow transition to turbulence occurs in the separated shear 
layer. The shear layer gains momentum from the free-stream and reattaches to the 
wall. Separation bubble length is inversely proportional to FSTI. 


3 Transition model coupling 

The V-model is not able to compute turbulence. Instead it determines the transi¬ 
tion threshold region. For this reason and as previously mentioned, the V-model 
transition closure was coupled to the Spalart-Allmaras turbulence model. Tran¬ 
sition onset prediction is performed by computing the viscosity induced by the 
predicted pre-transitional u v' values described throughout this work. The modus- 
operandi of the V-SA model is depicted in Figj2] 


4 Mechanical model approximation rational 


The rational behind the development of the transition V-model is herein pre¬ 
sented, including the flow physics on which it is supported. Before the mechanical 
approximation disclosure, some considerations need to be taken into account and 
explained. It is here assumed that pre-transitional turbulence is isotropic in a 
strain-less free-stream, in the sense that k x = k y = k z = k p . This is in agree¬ 
ment with the flow physics of transition such as presented in the work of |Mayle| 
and Schulz jl99 7|. However, under the effect of flow shear the model will predict 


small pre-transitional negative values of u v related to non-isotropic turbulence 
conditions. A bi-dimensional analysis is here described. Nevertheless the V-model 
transition closure is applicable to three-dimensional cases as will be later pre¬ 
sented. This is the case since the main three orthogonal shear deformation planes 
are all accounted for in the computation of the local shear magnitude. As such, a 
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global effect of three-dimensionality is taken into account in the computation of 
the pre-transitional uv . 


4.1 Pre-transitional turbulent kinetic energy considerations 


Pre-transition velocity fluctuations wave forms, for a specific frequency, have sel¬ 
dom a regular shape. Although this is true for most cases, the modeling of the 
pre-transition region, see works such as Jacobs and Durbin [1998 , requires both 
discrete and continuum modes. Citing Jacobs and Durbin [l 998| , ’’The eigenso- 
lutions to the Orr-Sommerfeld equation in an unbounded domain are classified 
into two spectra: the first is a finite set of discrete modes; the second is an in¬ 
finite continuum of modes. The latter are weakly damped and are irrelevant to 
classical linear stability analysis. Unstable modes are only members of the dis¬ 
crete spectrum. “ As previously stated, the modes used in classical linear stability 
analysis are the discrete spectrum components. However, the present transition 
model will attempt to model the effects on the continuum spectrum of modes. 
Citing Jacobs and Durbin [1998 , ’’The eigenfunction of the discrete modes de¬ 


cays exponentially with distance above the boundary layer. The eigenfunction of 
the continuous modes is sinusoidal in that region.” Therefore and in order to 
simplify the following exposure a sinusoidal wave shape was considered to model 
the free-stream pre-transition continuum spectrum. It is here assumed that si¬ 
nusoidal wave forms represent the time evolution of velocity fluctuations due to 
pre-transitional turbulence. Admitting that a particle is stuck inside one of these 
special pre-transitional turbulent vortices, its movement follows that of the vortex. 
Considering then a cross sectional plane of the bi-dimensional vortex, the equa¬ 
tions of motion for the particle imprisoned in the small pre-transitional vortex can 


be obtained. The equations of motion for x and y are then defined by ([lj) and ([2| 
respectively. These were obtained considering as a frame of reference the center of 
the pre-transitional vortex itself. 


X — Tavg COS 


V = -Tavg COS 



(1) 

( 2 ) 


The time derivative of the latter wil 
tuations u and v represented by ([i 


introduce the equations of the velocity fluc- 


31 and |4| 


u 


v 


= Tavgto Sm 

= r avg OJ Sin 



As already mentioned, it should be noted that these last four laws of motion are 
deduced assuming a frame of reference of the vortex itself. From the presented 
equations ([l]|4]), consideration of the a values must be taken. It can be seen that 
for a = ^ these laws of motion describe a circular motion as shown in Fig(3j 
This motion can be interpreted as a circular non-deformed pre-transitional vortex. 
Nonetheless, for a = f and a = ^ the described motion will be elliptical. The 













Transition V-model 


9 



Fig. 3 Particle trajectory describing a circular path through time. 


limiting cases are obtained for a = 0 and a = n, where although the motion is 
periodical, it is also linear. 

The pre-transitional vortex will have its turbulent kinetic energy. This can 
be related to a more simple definition of kinetic energy. In order to perform the 
analogy, we take another look at the clinging particle under the effect of the circular 
pre-transitional vortex. This circular motion analogy is considered under a bi- 
dimensional plane coincident with the xy Cartesian frame of reference positioned 
within the vortex center. Considering Fig|4j the animated particle in point P will 
have a certain amount of specific kinetic energy defined as, e Cp = U 2 /2. 

This same particle will have a mean turbulent kinetic energy since it follows 
the pre-transitional vortex rotational motion. The latter is defined as the sum 
of the mean turbulent kinetic energy along xx’s axis and the mean turbulent 
kinetic energy along yy’s axis. Assuming a circular vortex, and taking this into 
consideration, through § , a relation can be obtained in order to characterize the 
value of rotational velocity and radius of the pre-transitional vortex. The resulting 
relation is then <§. 


kx T" ky 


1-r-r 1 -r-r 

-u u + -V V = e Cp 


(iV avg ) 


( 5 ) 


\ 2 (kx T — Ldz'Tavg • (6) 

In <§> ojz represents the vortex rotational speed and r avg is the average pre- 
transitional vortex radius in the xy cross section plane. 

An equivalence of turbulent kinetic energy and kinetic energy was just per¬ 
formed. Besides the assumption of a circular pre-transitional vortex, up until this 
point no approximation has been used. However, in order to close this particular 
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Fig. 4 The pre-transitional turbulent vortices can be generated at the interface between 
the free-s tream turbulence and the pre-transitional laminar boundary layer due to ” inactive 
motion |Bradshaw| 1 994| . It is here presented a cut-section of a circular pre-transitional 
vortex. 


system an approximation must be performed. Either the average radius of the vor¬ 
tices is approximated or the rotational velocity of the vortices, oj z . The approach 
of rotational velocity approximation was chosen. The used relation was defined on 
the assumption that, rotational velocity of the pre-transitional turbulent vortices 
should be proportional to its pre-transitional turbulent kinetic energy. Also, the 
rotational velocity of pre-transitional vortices should be inversely proportional to 
the flow kinematic viscosity. This is further related to the fact that for a fixed 
turbulent large scale, the small turbulence sizes, such as the Kolmogorov scales, 
are reduced with increasing Reynolds number flows. The latter is also related to an 
increase of the turbulent vortices rotational velocities due to angular momentum 
conservation. Therefore, the lower the fluid kinematic viscosity, the higher the flow 
Reynolds number and consequently the higher the turbulent vortices rotational ve¬ 
locities. Besides this, dimensional analysis of the used relation confirms its correct 
dimensional characteristics. The selected approximation is then presented in 0 - 


Wz 


kx + ky 


( 7 ) 


Throughout the remaining model exposure ui z = w. Therefore, since k x + k y = 2k p , 
the value of to is now computed using the relation presented in ([8|. 


w = 


2 k p 
v 


( 8 ) 


From the latter relation and equation |6|, the value of r aV g is calculated in 


0 . 


T avg — 



( 9 ) 


In order to compute terms such as u u or u v' a time average of these fluctuating 
values must be performed according to equation (101. 
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= lim — 
T —>-oo T 


u u dt. 


( 10 ) 


It was then assumed a sinusoidal function for these fluctuating velocities. This 
purports that these fluctuations have a periodicity. The latter implies that a finite 
value for T is possible. A plausible value for T could be the periodicity of these 


sinusoidal functions. The latter assumption results in (111. 


u u = — / u u dt. 
27T L 


( 11 ) 


Using the explicit formulations of u and v , the calculation of u u and u v 


should be done according to (12) and (13 I respectively. 


2tt 


ravgOJ sin 


in (u*+f) 


dt. 


( 12 ) 


U V = 


27r 


r avg Lo sin (ut + ^ r aV gW sin (ut - ^ dt. (13) 


The u u value is directly obtained through (14). It can be seen that for 


isotropic pre-transitional turbulence, k x = k y , equation (141 is in accordance with 
the kinetic energy relations of ©> and ©. 


rLgU 2 


(14) 


The u v value is then calculated using relation (15). It should be noted that 


Tavg is the mean radius of the undeformed pre-transitional vortex. As can be seen 
in (15), there is a dependence with a. This a represents the phase shift between 


the two velocity components u and v . It also represents the deformation angle of 
the pre-transitional vortex. For a equal to ^rad or 90° the pre-transitional vortex 
has a circular undeformed shape. Also, for this value of a , according to (151, 


the Reynolds shear stresses, u v ', will be zero for an undeformed pre-transitional 
vortex. 


u v' 


rl vg u> 2 cos (a) 
2 


(15) 


4.2 Mechanical model for pre-transitional turbulent velocity fluctuation 
components under mean shear 

The fact that u v has a trend to present negative values under shear influence is 
commonly accepted. The reasoning begins by considering a no-slip wall constrained 
velocity profile. Under this scenario, a positive vertical velocity fluctuation away 
from the wall, such that v > 0, will induce a reduction of the flow momentum. 
The excited particle tends to maintain its momentum, this is lower than the new 
surrounding fluid due to the presence of the wall no-slip condition. This will then 
imply that u < 0, since the relative velocity of the new low momentum fluid is 
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Fig. 5 The deformation of a large scale vortical structure due to mean flow shear assumes the 
presented shape. Also, this is the expected pre-transitional vortex deformation shape under 
mean flow shear. 


lower than the surroundings. The reverse case is also applicable, that is for v < 0 
a value of u > 0 will follow. Therefore u v has a negative value. 

Consider then the case where an initially circular pre-transitional vortex is 
deformed by mean shear as in Fig[5] It must be noted that the presented schematic 
only demonstrates what is to be expected, not what the mechanical model predicts. 
Due to shear deformation, there will be an alteration of path curvature along the 
vortex surface. The centrifugal force distribution along the pre-transitional vortex 


will change accordingly. These centrifugal forces defined in (161, will act as pseudo 


non-linear springs of the pre-transitional vortex. These will act on the vortex when 
shear is present. The effect of the vortex deformation is reflected on the centrifugal 
forces computation, (16), through the local curvature radius, ri oca i- 


F'cent — 


U 2 pV 

'local 


(16) 


In order to compute the centrifugal forces, (16), the local radius of curvature, 
riocai , needs to be calculated. The change of the vortex curvature is dependent on 
the predicted deformation angle, a. As such, a relation between the local radius 
of curvature and the angle of deformation of the system is required. The latter 
relation was developed so as to deliver the requirements in, 


ric 


i.l oo : 


(17) 


& — 0 = ^ > 'f'local — Oj 


Q — 2 = * > ’ ^local — Tavg • 


(18) 

(19) 


The requirement presented in (171 is representative of an absolute flat vortex 
with orientation depicted in Fig[9| The following requirement in ( |18[ ) represents 
the scenario of a perfectly flat vortex with orientation depicted in Fig |5j The 
final requirement in (191 represents the undeformed circular vortex. Also, the 


calculation of local radius will depend on the average radius of the vortex and 
its deformation angle, a. The developed relation was then, 


Tlocal — Tavg~ 


(20) 


This relation enforces the requirements in (17p9|. 
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The shear force acting upon the vortex is defined in (21). In the latter, S, is 


mean flow shear magnitude. The A and V terms in both the equations (161 and 


(211, represent surface area and volume of the pre-transitional vortex respectively. 


These are defined in (22 I and (231. The l is here considered as the length of the 


pre-transitional vortex. A value for this is not required since l will cancel itself in 
the upcoming model development. 


ear — (M H - A, 

(21) 

A — < 2i'KTavgl i 

(22) 

V 'Kf'avgl' 

(23) 


This non-linear spring feature of the centrifugal forces can be modeled by a 
relatively simple mechanical system shown in Fig|6j The non-linear spring mechan¬ 
ical behavior is given by the centrifugal force physical characteristics as disclosed 
in (|16|). The ri oca i represents the local curvature radius of the pre-transitional 


vortex. In Figj6j the left picture represents the expected statistical mean shape of 
a pre-transitional vortex under the effect of mean flow shear. Normally the shear 
tensor main axes are aligned 45 ° with the local flow direction. The presented a 
and (3 angles in this left picture of Fig[6] are representative of the vortex defor¬ 
mation angles located in the directions of the shear stress tensor major axis and 
minor axis. Also, this a angle represents the phase shift in the motion equations 
( l][4). The right schematic in Fig j6] is the approximation of the continuum vortex 
shown in the left picture by a discrete system composed of four elements. In this 
schematic, the a and /3 angles are the half values of the original a and /? angles. 


Therefore, the angle relations in this problem are shown in (24) and (25). 


a = 2a 


(24) 


=- a 


(25) 


The solution of the mechanical dynamic problem presented in the right schematic 
of Fig (6] can be simplified by considering only one quarter of the system. This 
is possible due to the double symmetry along the shear stress tensor major axis 
and minor axis coincident with the axes of the ellipse resulting from the deformed 
circular pre-transitional vortex. With this in mind, the final mechanical model 
approximation is then shown in Fig[7] The presented orientation is in accordance 
to Fig(5j 

The attempt to solve the mechanical system in Fig[7]will produce a relation 
between vortex deformation and mean shear. In order to solve the mechanical 
problem presented in Fig[7J classical mechanical system solving procedures must 
be taken. First an equilibrium of the system moments in relation to point ”C” is 


in order. The resulting equation is presented in (261. 


SM C = R.A-^cos(a ) -I- F centA -sin(a ) - R B -sin(a ) - F centB ^cos(a ) = 0. 

(26) 
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B 



Fig. 6 The pre-transitional boundary layer oscillati ons due to shear known as KlebanofF modes 
were first observed and named ’’breathing modes“ by|KlebanofF 1971]. This behavior of the pre- 
transitional turbulent vortices can be accounted for by the mechanical model approximation. 
The depicted mechanical model approximation makes use of a fictitious non-linear spring 
analogy to describe the internal forces acting on the vortex. 



Fig. 7 The shear force, S, effect on pre-transitional eddies is approximated by this mechanical 
model for one quarter of a circular vortex. The mechanical approximation depicted here is 
applied to pre-transitional turbulent vortices in the upper zones of the laminar boundary 
layer. 


The shear force depicted in FigjT] by an S, is decomposed in the x and y axes 
directions resulting in S x and S y presented in (27) and (28) respectively. 


S x = S sin(a ). 


(27) 


S y = Scos(a ). (28) 

From this, the following required equation is the equilibrium of forces in the y axis 


direction. The obtained equation is disclosed in (29). 


SF y = Rb + F cent A + S y = 0 Rb = —FcentA — S y . 


(29) 
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The final required system equation is the x axis direction forces equilibrium equa¬ 
tion. The resulting equation is shown in (301. 

FF X = —Ra — FcentB + Sx = 0 Ra = — F cerl tB + S x . 


(30) 


equation (26) will result in (31). 


Substitution of the obtained relations (|29|) and (|30|) in the moment equilibrium 

h 


h 


h 


h 


h 


h 


-F cent ByCos(a )+S x —cos(a )+F centA -sin(a )+F centA -sin(a )+S y -sin(a )-F centB -cos(a ) = 0 

(31) 

(32) 


From the relations (|27j) and (281, equation (31) turns to (32). 

F cent Ahsin(a ) - F cent Bhcos(a ) + S hsin(a )cos(a ) = 0. 
From this, the following step is presented in, 

F cer itBCOs((X ) F 'cent ASin(a ) 


sin(a')cos(a') 


= S. 


(33) 


“A” and “B” are defined in, 

FcentA = 


According to equations (161, (20), (24) and (25), the centrifugal forces in points 


1 U 2 pV tt — 2a 

4 Tavg 2q 

(34) 

1 U 2 pV 2 a 

4 'f'avg 7T 2a 

(35) 


The value of the shear force depicted by an S in Fig[T] is equal to a quarter of the 
total shear force acting upon the pre-transitional vortex defined in (21). Therefore, 
S is defined in, 

S = \(p + p T )SA, (36) 

As such, from these relations the resulting equation from (33 I is disclose in, 

(rfiv)»“(“)-( I iv !: )™(“')' 


U z pV 


sin (a') cos (a') 


= (p + p,T)SA. (37) 


From the mean linear velocity of the pre-transitional vortex defined in (381 and 


from the definition of area, “A“, and volume, “V“, defined in (22) and (23) re¬ 


spectively, the ratio between pre-transitional vortex acting shear and centrifugal 


forces is given in (39). 


U — OJTavgi 

(38) 

, • F s hear {y H” )2*S 

ratio = — - = --—. 

(39) 

recent (yJT'avg') 


The solution of this system in terms of angle a is then, 

(jig) c ° s 0*') - (gf£) sin ( a ) 

sin ( a ') cos (a') 


= ratio. 


(40) 
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The evolution of ratio to a of the exact solution (40) is presented in Figjs} 
According to relation (401 and observing Fig[8j an unexpected evolution of the 
angle a is predicted by the mechanical model. Based on ( |4o[ ), this solution of the 
mechanical model approximation used in this work predicted an increase of the 
angle a with increasing shear. This in turn means that according to (151 and (241, 
the calculated values of u v will be negative. The a angle increases from J to its 
limiting value 7r. This also signifies that the correct pre-transitional vortex shape 
under the effect of shear is the one presented in Fig [9] Since the relation (40) is 
far too complex to be use in a numerical transition closure, a simpler function is 
used instead. The used function is a good approximation of the exact relation and 
is presented in, 


tan 


1.8 


(“ 



= ratio. 


(41) 


The used mechanical approximation function in the transition V-model numer¬ 
ical implementation is here presented in (42 I. This is done in order to simplify 
implementation attempts by the reader. 


a = mm 


/ .. \ 1.42857 

—i / ratio \ 7T 7T 

tan ( - -I—, — 

1.8 J 4’ 2 


(42) 


In (42), the minimum function is applied since the maximum value of a is 
This correspond to a perfectly flat pre-transitional vortex deformation angle. For 
values of ratio = 0 the minimum value of a is J. This represents an undeformed 
circular pre-transitional vortex. 


















Transition V-model 


17 


S„ 




Fig. 9 The derived relation ( |40| ) for co mpu tation of the a deformation angle is used to 
calculate the a deformation angle through ( |24| ). The former relation predicts an a angle increase 
with shear, instead of the expected angle decrease a s de picted in Fig[5] As the mean flow shear 
increases so does the ratio calculated according to ( |39| ). The predicted pre-transitional vortex 
deformation under mean flow shear is presented here. 


5 The transition V-model transport equation 


The Reynolds stress transport equations were first closed by Rotta [195 1 , setting 
the foundations for second order turbulence models. The transport equations are 
here presented in (43) as presented by Wilcox 1994 . 


dT ij ,rr dT ij 

dt k dxk 


— T~ik 


dU±_ 

dxi. 


T ;ik 


dUi 

dx k 


+ 7 Ui j + 


dx k 


dr z . 
/-- 

dxi 


+ Ci 


ijk 


, (43) 


M dx k dx k 


(44) 


, / du- du 

n ij =p hr + ir 1 


dxj 


dx-i 


(45) 


Cijk = pu^u'-Uf. +p'u' i 5 jk +p'u'Afc- 


(46) 


In (431, Tij = —pu'iUj, represents the Reynolds stress tensor. The presented 
terms on the right hand side of (43), are from left to right, two terms of production, 


one for turbulent dissipation presented in (441, another for pressure-strain shown in 


45k and finally the molecular and turbulent diffusion. This last term is presented 


in (|46j). Considering the normal Reynolds stresses transport equations, that is, the 
trace components of the stress tensor, the number of terms is reduced resulting 


in (47). Through the condition imposed by the equation of continuity applied to 


turbulent velocity fluctuations the pressure-strain term is null. 


dm dm _ „ dU l du i du i d 

—h Uk y;— — ~2 T ik ,,-b — ,,-b y,— 

at oxk oxk oxk oxk oxk 


dm . — i—r-r -— x 

+ P U i U i U k + 2 P U i°ik 


_ (47) 

Since ki = considering then the relation that, m = — 2 pfcj, we obtain 

the following transport equation for turbulent kinetic energy in (481. 
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dkj dkj 

dt k dxk 


— dUi 
lk dx k 


^du\ du'i d 

dx k dx k dx k 


dki 

/ - 

dx k 


-pu'^u'k -p'u'iSik 


(48) 

The proposed transition model pre-transition turbulent kinetic energy trans¬ 


port equation components are based on some terms of equation (481. The tran¬ 


sition closure production term was obtained by analyzing the first term in (481. 


This analysis will be performed considering a flat-plate flow far from its leading 
edge region. Within the shear region of the flow, or, u'v', will be negative. 

The shear value of dU\ /dx 2 will be positive. Thus this term is a production term 
in mean flow shear conditions. The mean flow shear magnitude, S, is used as the 
mean flow property that drives pre-transitional turbulence production. It must be 
noted that for clarity reasons, a term was omitted in the previously presented u v' 


calculation in (151. The term is a function responsible for the evaluation of the 


mean flow and turbulent scales proximity, Fs- This is defined in (49). 


Fs = 1 — min ( mag 


( Ssize T a 


\max ( S. 


sxze.’t • avg 


The used mean flow shear scale function, S s ize, is defined in (50). 


Ssize — 

The u'v' computation is then performed according to, 

—!—r (ravgU)) 2 cos (a) „ 

U V = --— --— -Fa. 

2 

The production term used is presented in, 


Prodk = —C Pk Su'v'. 


(49) 


(50) 


(51) 


(52) 


In the latter equation, Cp k is the V-model calibration constant equal to 0.835. 


The second term on the RHS of (48) will be used to obtain the dissipation 


term of the presented transition model transport equation. Calculation of the 
cross partial derivatives from the latter term reveals a fundamental result. In 
order to simplify the partial derivatives calculation, the following assumption is 
required. The equations (JT]) and ([ 2 ]) definitions impose a phase shift between x 
and y . For an undeformed pre-transitional vortex, implying a = ^, and summing 
a phase shift value of J to both equations, (jlj) and (|2j), we obtain the same phase 
shift imposition through (55) and (56) respectively. A detailed explanation of the 


process is presented in (53) and (541. 


x = —r a 


V = ~r a 


, cos 


Ujt+ 2 + l 


, sin (wt), 


ut - + — J = —r avg cos (ut ). 

= Vavg sin (cot) , 


(53) 

(54) 

(55) 
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y = ravg cos (cot) . 


(56) 


The new velocity fluctuations u and v are then represented by (57) and (581 
respectively. 


u = r avg co cos (cut ), 


(57) 


v = r avg co sin (cot). 


(58) 


Considering then a cross partial derivative where the length scale of dy is compa¬ 
rable to the length scale of dy resulting in, 


(59) 


du du dr avg co cos (cot) 
dy dy' dy' 

It can be observed that apparently the latter partial derivative can not be easily 
solved. However, the cot dependence can be written as, 


y = —ravg cos (cot) 44- cot = arcos 


By doing so the cross partial derivative becomes then, 


du 

dy' 


dvavg COS 


( arcos (fd;)) dcor a 


dy' 


dy' 


The resulting cross partial derivative is presented in (62 I. 


du 

dy' 


dcoy 

dy' 


(60) 


(61) 


(62) 


Therefore the cross partial derivative of the dissipation term (441 is equal to a 


form of rotational velocity of the undeformed pre-transitional turbulent vortices, 
co. It must be noted that due to the assumption of dy ss dy , a scale proximity 
between turbulent and mean flow must exist. The non-cross partial derivatives of 


(441 were not considered. The resulting dissipation term in the V-model transport 


equation is then presented in (63). It must be noted that although the cross partial 
derivative, (621, of the dissipation term, (441, has a negative sign, the dissipation 


term makes use of the square of this cross partial derivative. 


Destk„ = vf2 Fn ~ v 


du- du ' 


dx k dx k ' 


(63) 


Similar to the used ’’shear-sheltering” effect function presented in the work of 
Walters and Cokljat [2008 , mean flow vorticity, 17, is used as a form of rotational 
velocity instead of the pre-transitional vortex rotational speed, co, for the destruc¬ 
tion effect within the boundary layer. Citing from the work of[Walters and Cokljat] 
||008|, “Shear-sheltering refers to the damping of turbulence dynamics that occurs 
in thin regions of high vorticity...”. This is also done in order to allow the mean 
flow characteristics to control transition onset as depicted in Fig[2] Also, mean flow 
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vorticity can be related to stabilization effects of turbulence, which can be inter¬ 
preted as a turbulence sink. Such can be observed in relaminarization experiments 
of turbulent flow inside tube coils as presented by the work of Viswan ath et al.] 
[1978] . This was also calculated with DNS, Noorani et al. [2013 . The Fq function 
assures the aforementioned assumption of dy ss dy . Fq can also be interpreted as 


a measure of the mean and turbulent flow scale proximity. This is defined in (641. 


Fq = 1 — min mag 


( f^size T a 


,1 


\^maX size: Tavg 

The used mean flow vorticity scale function, F2 s i ze , is defined in (651. 


FI size. — \ ^ 


(64) 


(65) 


The turbulent diffusion component of the transport equation was chosen to be 
equal to a common turbulent kinetic energy transport equation such as the model 
of Craft et al. 1996 . The resulting pre-transitional turbulent kinetic energy trans¬ 
port equation is then presented in, 


Dk p 

~Dt 


= Prodb - DesU + 


d 

dxj 


( i \ dk p 

( " + " T >57 


( 66 ) 


From the transition V-model, the pre-transition turbulent kinetic energy will apply 


an induced viscosity defined in (67). 


^Tuv — 


— U V 


mag (Vf7) 


(67) 


This will then be the small effect of the pre-transitional turbulent kinetic energy 
on the viscosity within the pre-transitional region of the laminar boundary layer. 


6 Transition V-model and Spalart-Allmaras turbulence model coupling 

The Spalart-Allmaras is a well known one-equation turbulence model. It was first 
presented by Spalart and Allmaras 1994 . As stated before, the transition V-model 
was coupled to the Spalart-Allmaras turbulence model. This was performed by 
simply adding a control function in the production term of the turbulence model. 
This control function depends upon the ratio value obtained by the mechanical 
model approximation in (391. The performed implementation is shown in (681. 


Dvt 

Dt 


= CbiSFrexp (—Cv~ T ratio) — Cwij 


/ 

\\ — 

( - , v\ 

dv T 


dxj 

[v T F-) 

dxj 


( 68 ) 


As can be seen in the first component of the RHS of (681, the control function is 


an exponential term. The C„~ T constant is equal to 0.5. Also in order to account for 
the effect of the pre-transitional viscosity, (671, the Spalart-Allmaras turbulence 


model production term is slightly changed to, 

Fl 


S = 


+ 


Uj 1 


Vi M ) 2 


(69) 



































Transition V-model 


21 


Table 1 ERCOFTAC ZPG Flat-Plate Upstream Condition 


Case 

Tu {%) 

U {m/s) 

T3A 

3.0 

5.4 


For brevity’s sake, further description of the Spalart-Allmaras model is omitted 
here. The remaining Spalart-Allmaras model is briefly described inlRumsey et al.l 
|001 . 


The total turbulent viscosity that the V-SA transition model predicts is a sum 
of the obtained turbulent viscosity from the SA model and the transition V-model. 
This is presented in (701. 


VT = l'Tuv + fvlVT- 


(70) 


7 Results and Discussion 


All of the obtained results for the V-SA and SA models were calculated using the 
open-source software OpenFoam. The considered cases were all steady-state and 
incompressible. These were performed with a constant density of p = 1.2 (fc<?/m 3 ). 
The results calculated with OpenFOAM were run with a pressure based solver 
SIMPLE, linear discretization for laplacian terms and LUST discretization scheme 
for all possible divergence terms. In OpenFOAM the divergence term, (711, requires 
a linear discretization. 


div((nuEf / * dev(grad(U).T ()))). (71) 

All of the obtained results with the empirical transition model 7 — R e g were 
computed using Ansys Fluent 13.0. The used discretization scheme was the second 
order upwind, SOU, and for pressure the linear setting was applied. 

The V-SA wall boundary conditions for all tested cases were Dirichlet condi¬ 
tions for pre-transition turbulent kinetic energy, that is, k p = 0 at the wall. 

As stated in the Introduction section, the first comparison is performed by using 
the results from ERCOFTAC obtained by the experimental work of jCoupland] 
[199 Ob I, for one of the zero-pressure-gradient, ZPG, flat-plate test cases. The tested 
case is then the T3A. The upstream conditions for this test case is presented in 
table |T| 

For the zero pressure gradient test case the flat-plate mesh used was structured 
and had y + values below 0.1. The flat-plate had 1.7 meters of extension with 200 
mesh points over its surface. These clustered near the leading edge of the plate. 
The leading edge had a curvature radius of 0.002 meters. Along the leading edge 
the mesh had 30 mesh points. The wall perpendicular spacing of the first layer of 
cells over the flat-plate was 1 x 10 -5 meters. The velocity inlet was located at 0.15 
meters from the leading edge. This short extension had 110 mesh nodes. The top 
surface was located at 0.15 meters above the flat-plate. This vertical length had 
110 mesh points along it. 

The inlet boundary conditions for this flat-plate ZPG ERCOFTAC test case 
is presented in table [ 2 ] 

The results for the flat-plate T3A ERCOFTAC test case are presented in 
Fig|10| These include a mesh refinement validation. This was performed using an 
equal mesh with double number of nodes over the entire computational domain. 
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Table 2 ERCOFTAC ZPG Flat-Plate Inlet Boundary Conditions 


Case 

kp(m z / s z ) 

U(m/s) 

VT /v 

T3A 

0.039366 

5.4 

3 



Fig. 10 Comparison of experimental data from ERCOFTAC T3A flat-plate test case skin- 
friction coefficient distribution with the non-transition SA closure and the transition V-SA 
model. The used structured computational grids for the mesh refinement study had the re¬ 
spective sizes of 110 X 340 and 220 X 680 nodes. 


As can be observed, although the Spalart-Allmaras model is able to calculate 
a certain transition, it predicts transition onset too early. The V-SA transition 
model improves on the transition onset prediction, demonstrating a more correct 
behavior. The V-SA mesh refinement study indicates that the results are mesh in¬ 
dependent. As already mentioned, the V-model transition closure calculates small 
pre-transitional negative values of u v . This is shown for the region near the lead¬ 
ing edge and for the transition onset zone of the T3A test case in FigflTl As can be 
seen, when the mean flow characteristics are ideal the V-model predicts the pierc¬ 
ing of the laminar boundary layer by the pre-transitional u v . This then activates 
the SA turbulence model production term. 

The proposed mechanical model approximation for pre-transitional vortex de¬ 
formation predicts the distribution of u v in the pre-transitional boundary layer 
region. A comparison is performed between the experimental ERCOFTAC database 
of u v and the V-model predicted values. These comparisons are performed from 
the leading edge to the transition onset point. The last comparison is one station 
in the middle of the transition region. The results are presented in Figsjrg [T3J 
|14[ |15| |l6| and [IT] These represent the axial positions in meters corresponding to 
0.095, 0.195, 0.295, 0.395, 0.495 and the transition section 0.595 respectively. The 
latter axial positions correspond to the Reynolds numbers of 3.24 x 10 4 , 6.70 x 10 4 , 
10.06 x 10 4 , 13.48 x 10 4 , 16.92 x 10 4 and the transition section 20.35 x 10 4 respec¬ 
tively. 

As can be seen, the V-model predicts an overshoot of u v near the leading 
edge. However, the results greatly improve along the pre-transitional region. The 
results in Fig|17| corresponding to a station in the middle of the transition length 
region, show that the model predictions deviate from the measured values. 
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0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 


Fig. 11 Contour map of u v values with detailed results near the leading edge and transition 
onset zone of the T3A flat-plate test case. 



Fig. 12 Comparison of ERCOFTAC flat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.095 meters or Rex of 3.24 x 10 4 . The 
ne gative pre-transit ional u v values were related to ” splat-mechanism “ or ” inactive-motion “ 
by |Bradshaw| [1994] . 
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Fig. 13 Comparison of ERCOFTAC flat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.195 meters or Rex of 6.70 x 10 4 . 



Fig. 14 Comparison of ERCOFTAC flat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.295 meters or Rex of 10.06 x 10 4 . 


Table 3 ERCOFTAC Pressure Gradient Flat-Plate Upstream Conditions 


Case 

Tu(%) 

U(m/s) 

T3C3 

3.0 

3.7 


For the T3C test case the ERCOFTAC flat-plate pressure gradient transition 
experimental data from Coupland 1 1990a| was used. The tested case was the T3C3. 
The upstream conditions for this test case are presented in table [3] The experimen¬ 
tal pressure gradient flat-plate test case was performed with a structured mesh. 
The bottom part of it was equal to the ZPG T3A test case mesh. The upper 
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Fig. 15 Comparison of ERCOFTAC fiat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.395 meters or Rex of 13.48 X 10 4 . 



Fig. 16 Comparison of ERCOFTAC fiat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.495 meters or Rex of 16.92 X 10 4 . 


component of the mesh was a curved wall surface designed to obtained the exper¬ 
imentally measured free-stream velocity variations. The curved surface had 200 
mesh points along it. The mesh points cluster near this surface and the spacing of 
the first layer of cells was 1 x 10 -4 meters. This top section of the mesh had 30 
nodes along its vertical connectors. Summing this with the previous bottom part 
of the mesh makes a total of 140 mesh points in the vertical direction over the 
flat-plate. The inlet of the complete mesh had a vertical length of 0.3 meters. The 
inlet boundary conditions for the T3C3 ERCOFTAC test case are presented in 
table [4] The results for the flat-plate T3C3 ERCOFTAC test case are presented 
in FiglU As can be observed, the V-SA model’s transition onset prediction is 
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Fig. 17 Comparison of ERCOFTAC flat-plate T3A experimental u v values with those pre¬ 
dicted by the transition V-model in the axial position of 0.595 meters or Rex of 20.35 X 10 4 . 


Table 4 ERCOFTAC Pressure Gradient Flat-Plate Inlet Boundary Conditions 


Case 

k p (m 2 / s 2 ) 

U(m/s) 

i/t/v 

T3C3 

0.01848 

3.7 

3 



Fig. 18 Comparison of experimental data from ERCOFTAC T3C3 flat-plate test case skin- 
friction coefficient distribution with the non-transition SA closure and the transition V-SA 
model. 


in accordance with the experimental data. The SA turbulence model predicts an 
early transition onset. 


As mentioned earlier, the ERCOFTAC T3L test cases of separation induced 
transition were considered for validation purposes. The importance of such vali- 


dation is made clear in the work of 

dadzic and Hanjalic 

11999 . The experimental 

results from 

Coupland and Brierley 

11996 , were used for the present validation. 


The used computational mesh was structured and had y + values below 0.1 along 
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Table 5 ERCOFTAC T3L Flat-Plate Upstream Conditions 


Case 

Tu(%) 

U (m/s) 

T3L1 

0.2 

5.0 

T3L3 

2.3 

5.0 

T3L5 

2.3 

2.5 


Table 6 ERCOFTAC T3L Flat-Plate Inlet Boundary Conditions 


Case 

k p (m 2 / s z ) 

U(m/s) 

l/J? / V 

T3L1 

0.00015 

5.0 

5 

T3L3 

0.019838 

5.0 

5 

T3L5 

0.004959 

2.5 

5 


the whole flat-plate surface. The flat-plate had 1.7 meters of extension with 200 
mesh points over its surface. These clustered near the leading edge of the plate. The 
leading edge had a curvature radius of 0.005 meters, which is in accordance with 
the experimental setup. Along the leading edge the mesh had 60 mesh points. 
The wall perpendicular spacing of the first layer of cells over the flat-plate was 
1 x 10 -5 meters. The velocity inlet was located at 0.15 meters from the leading 
edge. This short extension had 110 mesh nodes. The top surface was located at 
0.15 meters above the flat-plate. This vertical length had 110 mesh points along 
it. The upstream conditions for the considered test cases are presented in table 
[5] The inlet boundary conditions for the flat-plate ERCOFTAC T3L test cases 
are presented in table [ 6 ] For these cases, separation induced transition will be 
analyzed. As such, the fluid kinematic viscosity used was the experimental value 
of v = 1.6 x 10 -5 (m 2 /s). As shown in table 5 the T3L1 test case represents 
the lower bound of turbulence intensity. The results for this case are presented in 
Fig|19| The flat-plate leading edge oscillations of skin-friction coefficient are due 
to laminar-boundary layer separation. The results show that the V-SA transition 
model has difficulty predicting transition onset under these turbulence and flow 
conditions. However, transition still occurs later on over the flat-plate. The SA 
model is able to handle these conditions better, determining flow transition at the 
correct location. 

The second separation induced transition test case considered has a higher tur¬ 
bulence intensity. The T3L3 ERCOFTAC flat-plate test case results are presented 
in Fig|20| As can be seen, the V-SA and SA models are able to determine the 
transition onset very close to the experimental data. Again the flat-plate leading 
edge oscillations of skin-friction coefficient are due to boundary layer separation. 
R should be noted that both tend to slightly delay transition onset prediction. 

Finally for the T3L5 test case the results presented in Fig]21| indicate that 
both the V-SA and SA closures are able to determine transition onset near the 
experimental measurements. For both cases the flat-plate leading edge oscillations 
of skin-friction coefficient are due to flow separation. 

The incompressible three-dimensional geometry of a 6:1 prolate-spheroid was 
computed using SA, V-SA and the empirical correlation transition model of Ansys 
Fluent, the 7 — R e g. These test cases were performed using the experimental data 
[1985 


of Kreplin et al. 


This was obtained through a personal communication 
with Dr. Kreplin. The latter experimental work has many test cases, however only 
some of these were considered for validation purposes. The selected three test cases 


















28 


Rui Vizinho et al. 



Fig. 19 Comparison of experimental data from ERCOFTAC T3L1 flat-plate test case skin- 
friction coefficient distribution with the non-transition SA closure and the transition V-SA 
model. The flat-plate leading edge oscillations of skin-friction coefficient are due to flow sepa¬ 
ration. 



Fig. 20 Comparison of experimental data from ERCOFTAC T3L3 flat-plate test case skin- 
friction coefficient distribution with the non-transition SA closure and the transition V-SA 
model. The flat-plate leading edge oscillations of skin-friction coefficient are due to flow sepa¬ 
ration. 


had significant transition effects. The upstream conditions for these test cases are 
presented in table [7] 

The experimental setup comprised of a 6:1 prolate-spheroid with a major and 
minor axis lengths of 2.4 and 0.4 meters respectively. For the considered test cases 
the experimental flow velocity was 45 (m/s). The experimental fluid kinematic vis¬ 
cosity was v = 1.7 x 10 -5 (m 2 /s). In order to perform validation using these test 
cases experimental data, an equal 3D geometry was used with the same dimen¬ 
sions. The used prolate-spheroid mesh was structured and had values below 
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Fig. 21 Comparison of experimental data from ERCOFTAC T3L5 flat-plate test case skin- 
friction coefficient distribution with the non-transition SA closure and the transition V-SA 
model. The flat-plate leading edge oscillations of skin-friction coefficient are due to flow sepa¬ 
ration. 


Table 7 6:1 Prolate Spheroid Upstream Conditions 


AoA 

Tu(%) 

Re 

5° 

0.1 

6.5 x 10 H 

15" 

0.1 

6.5 x 10 B 

30” 

0.1 

6.5 x 10 b 


Table 8 6:1 Prolate Spheroid Inlet Boundary Conditions 


AoA 

k p (rn z /s z ) 

U(m/s) 

V~t /v 

5" 

0.003038 

45 

3 

15" 

0.003038 

45 

3 

30° 

0.003038 

45 

3 


0.6 over the entire surface. Along the surface major axis the mesh had 400 com¬ 
putational nodes. In the azimuth orientation, the prolate-spheroid cross-section 
had 100 mesh nodes. The total value of grid points over the prolate-spheroid sur¬ 
face was 40000. The first layer of cells over the spheroid surface were distanced at 
1 x 10 -5 meters. The inlet boundary conditions for these test cases of transition 
under cross-flow effects over a 6:1 prolate-spheroid are presented in table | 8 j The 
applied inlet boundary conditions for the Fluent transition model 7 — R e g , were 
selected in order to simulate the same turbulence intensity in the free-stream as 
presented in table [7J The fluid kinematic viscosity used was the experimental value 
of v = 1.7 x 10 -5 (m 2 /s). 

Cut-section plots of skin-friction coefficient for the obtained results is per¬ 
formed. The cut-section plane is perpendicular to the 6:1 prolate-spheroid minor 
axis and contains the latter volume center point. This section cuts the 6:1 prolate- 
spheroid in the x-z plane as shown in Fig|22| 
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Fig. 22 X-Z cut-section of 6:1 prolate-spheroid for skin-friction coefficient plots. The cutting 
plane is perpendicular to the 6:1 prolate-spheroid minor axis and contains its origin point. The 
presented structured mesh has the size of 100 x 400 nodes. 



Kreplin et al, 


1985 


Fig. 23 Comparison of experimental data of 
contour map of 6:1 prolate-spheroid with AoA 5 “ with numerica 
SA closure and transition V-SA and 7 — R e e models. 


for skin-friction coefficient 
results of non-transitional 


The results for the 5° angle of attack, AoA, test case are presented in Fig|23[ 
The experimental data for the prolate-spheroid tips is not available. Thus these 
have been assigned with a zero value of skin-friction coefficient in all presented re¬ 
sults. However, it must be noted that for all AoA, the flow separates at the trailing 
edge of the prolate-spheroid. As can be seen, the SA closure determines transition 
onset right at the leading edge of the prolate-spheroid. The V-SA transition model 
predicts transition onset near to the experimental result. It should be noted that, 
the transition length of the V-SA is shorter than the experimental data. However 
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Fig. 24 Comparison of experimental skin-friction coefficient along a top X-Z cutting plane 
over the 6:1 prolate-spheroid with AoA 5"with numerical results of non-transitional SA closure 
and transition V-SA and 7 — R e g models. 


the transition onset position is quite close to the experimental measurements. The 
7 — R e o transition model predicts the transition onset correctly but the transition 
line along the surface has an incorrect angle. A skin-friction coefficient plot from a 
top x-z cutting plane parallel to the one presented in Fig[22]is presented in Fig|24| 
It is shown that the 7 — R e g empirical transition model is able to correctly predict 
transition onset as well as the late transition value of skin-friction coefficient. The 
central x-z cutting plane results are presented in Fig|25[ The presented results 
confirm the reliability of the V-SA model. 


The results for the AoA 15 ° test case shown in Fig|26| expose the fact that 
the SA model predicts transition at the leading edge of the geometry. The V-SA 
closure is able to predict transition onset near the experimental values although 
with some delay. Again the 7 — R e g transition model predicts the transition onset 
very close to the experimental data but the transition zone shape is incorrect. In 
the work of Seyfert and Krumbein| [2012 , similar numerical results were obtained 
with the 7 — R e g transition model for these last two test cases. The x-z cutting 
plane results are presented in Fig|27| As can be seen, the V-SA model predicts 
transition onset close to the experimental data. However, the predicted transition 
skin-friction coefficient peak value is lower than that of the experimental data. 
Also, for this test case a velocity profile normal to the prolate-spheroid surface was 
taken in order to observe the cross-flow effect. The velocity profile twist due to the 
present cross-flow conditions is shown in Fig|28| The evolution of flow streamlines 


and u v' iso-surfaces over the spheroid are presented in Fig 29 As can be observed, 


there is flow separation at the trailing edge of the spheroid. Also, the leading edge 
u v' iso-surfaces patterns are quite interesting. There seems to be two sets of u v' 
fluctuations over the spheroid nose. This can also be observed in the front view of 
the 6:1 prolate-spheroid u v' iso-surface patterns. These are presented in Fig 30 


In the fully turbulent flow region these patterns cease to exist. Instead a constant 
iso-surface covers the remaining extension of the prolate-spheroid. 
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x(m) 

Fig. 25 Comparison of experimental skin-friction coefficient along the X-Z cutting plane over 
the 6:1 prolate-spheroid with AoA 5“with numerical results of non-transitional SA closure and 
transition V-SA and 7 — R g q models. 


SA 


Experimental 


V-SA 






Cf 


D.QD75 

0.00715 

0.0068 

0.00645 

0.0061 

0.00575 

0.0054 

0.00505 

0.0047 

0.00435 

0.004 

0.00365 

0.0033 

0.00295 

0.0026 

0.00225 

0.0019 

0.00155 

0.0012 

0.00085 

0.0005 


Fig. 26 Comparison of experimental skin-friction coefficient contour map of 6:1 prolate- 
spheroid with AoA 15 ° with numerical results of non-transitional SA closure and transition 
V-SA and 7 — R g q models. In the regions marked by the letters A and B there is a severe 
cross-flow transition effect. 
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Fig. 27 Comparison of experimental skin-friction coefficient along the X-Z cutting plane over 
the 6:1 prolate-spheroid with AoA 15 ° with numerical results of non-transitional SA closure 
and transition V-SA and 7 — R e Q models. 


Y 



X 



Fig. 28 Cross-flow effect on velocity profile normal to the surface of the 6:1 prolate-spheroid 
with AoA 15 °. The V-SA transition model calculated skin-friction coefficient contour map is 
used as surface contour. 


The final 6:1 prolate-spheroid validation test case was performed with an AoA 
of 30 °. As can be seen in the results of Fig|3TJ the V-SA model’s transition on¬ 
set prediction is in accordance with the experimental measurements. The V-SA 
transition closure can even predict the saw-tooth shape behavior of the transition 
onset line very similar to the experimental data. The SA turbulence model predicts 
transition onset at the beginning of the prolate-spheroid. The 7 — R e g exhibits a 
similar behavior to the last two test cases presented here. The x-z cutting plane 
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Fig. 29 Side view of red colored transparent u v iso-surfaces with flow streamlines over 
the 6:1 prolate-spheroid with AoA 15 Top image represents the u v iso-surface equal to 
-0.02. Bottom image represents the u v iso-surface equal to -0.03. The V-SA transition model 
calculated skin-friction coefficient contour map is used as surface contour. 


results are presented in Fig|32( As shown the V-SA model is able to predict tran¬ 
sition onset near the experimental data and calculates the transition process with 
a saw-shape. The latter shape is due to strong cross-flow effects during turbulence 
transition. This effect is visible in the regions marked by the letters A and B in 

Fig|U 


8 Conclusions 


The rational of the transition V-model development was presented. A detailed 
description of the mechanical approximation model was performed. The transition 
model pre-transitional turbulent kinetic energy transport equation formulation 
was derived. 

The model was validated for eight benchmark test cases. The V-SA transition 
model was validated in the zero-pressure gradient flat-plate T3A test case. The V- 
model predicted values of u v in the pre-transition boundary layer were validated 
with experimental data from the flat-plate T3A ERCOFTAC database. As stated 


by Bradshaw 1994 , ’’the so-called inactive motion... is simple: the motion near the 
surface, ... results mainly from eddies actually generated near the surface, ... the 
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Fig. 30 Front view of red colored non-transparent u v iso-surfaces over the 6:1 prolate- 
spheroid with AoA 15 °. The image represents the u v iso-surface equal to -0.03. The V-SA 
transition model calculated skin-friction coefficient contour map is used as surface contour. 



Fig. 31 Comparison of experimental skin-friction coefficient contour map of 6:1 prolate- 
spheroid with AoA 30 ° with numerical results of non-transitional SA closure and transition 
V-SA and 7 — R g q models. In the regions marked by the letters A and B there is a severe 
cross-flow transition effect. 
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Fig. 32 Comparison of experimental skin-friction coefficient along the X-Z cutting plane over 
the 6:1 prolate-spheroid with AoA 30 ° with numerical results of non-transitional SA closure 
and transition V-SA and 7 — R e g models. 


contribution... to the shear stress -pu v' is small". The transition V-model predicts 
this small contribution. Also the V-SA model was validated with the flat-plate 
pressure-gradient ERCOFTAC T3C3 test case. The experimental ERCOFTAC 
T3L flat-plate test cases were used to further validate the V-SA model. These 
were used for validation of the V-model under separation induced transition. For 
the very low turbulence intensity case of Tu=0.2%, T3L1, the model predicted 
a delayed transition onset. This is possibly related to an excessive effect of the 
V-model destruction term. Since the free-stream turbulence intensity is very low, 
with low free-stream velocity, so is its turbulent kinetic energy. Therefore there 
seems to be an excessive sensibility to the destruction effect by the separation 
bubble vorticity. For all of the remaining T3L tested cases the model behaved 
correctly. 

The strengths of the V-SA model are verified on the transition under cross- 
flow effects test cases of the three-dimensional 6:1 prolate-spheroid geometry. It 
was observed that although the free-stream turbulence intensity for these cases is 
very low, Tu=0.1%, the model is able to predict transition onset patterns close 
to the experimental data. However, the transition length of the model is short in 
one of the tested cases. For the test case with a low angle of attack of 5 °, the 
V-SA model predicts a short transition length. Although the reason for the latter 
is unclear, it is suspected that an excessive pre-transition turbulent kinetic energy 
diffusion inside the boundary layer might be the reason for such short transition 
length. The rate of turbulence intermittency diffusion into the transition boundary 
layer has a major role determining transition length as shown in the work offDurbin 
[ 2012 [ . 

For further research the transition length control of the model should be in¬ 
vestigated as well as the delayed transition onset under the combination of very 
low free-stream turbulence intensity and low speed conditions. 
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